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ABSTRACT 

We calculate the theoretical evolution of the radii of all fourteen of the known transiting extrasolar 
giant planets (EGPs) for a variety of assumptions concerning atmospheric opacity, dense inner core 
masses, and possible internal power sources. We incorporate the effects of stellar irradiation and 
customize such effects for each EGP and star. Looking collectively at the family as a whole, we find 
that there are in fact two radius anomalies to be explained. Not only are the radii of a subset of 
the known transiting EGPs larger than expected from previous theory, but many of the other objects 
are smaller than the default theory would allow. We suggest that the larger EGPs can be explained 
by invoking enhanced atmospheric opacities that naturally retain internal heat. This explanation 
might obviate the necessity for an extra internal power source. We explain the smaller radii by the 
presence in perhaps all the known transiting EGPs of dense cores, such as have been inferred for 
Saturn and Jupiter. Importantly, we derive a rough correlation between the masses of our "best-fit" 
cores and the stellar metallicity that seems to buttress the core-accretion model of their formation. 
Though many caveats and uncertainties remain, the resulting comprehensive theory that incorporates 
enhanced-opacity atmospheres and dense cores is in reasonable accord with all the current structural 
data for the known transiting giant planets. 

Subject headings: stars: transits — planetary systems — planets and satellites: general 



1. INTRODUCTION 

Approximately two hundred extrasolar giant planets 
(EGPs) have to date been discovered by radial-velocity 
techniques'*. These data yield orbital properties and 
MpSin(i), where Mp and i are the planet mass and in- 
clination angle, respectively. However, for the subset of 
fourteen EGPs that are currently known to transit their 
primaries (Charbonneau et al. 2006a), the Mp— sin(i) 
degeneracy is broken and EGP radii (Rp) are measured 
as well. With Mp and Rp, an estimate of the (presum- 
ably) coeval stellar age, and a detailed theoretical model 
that includes the effects of stellar irradiation, the general 
theory of the structure, evolution, and atmospheres of ir- 
radiated close-in EGPs^ can be put to the test (Burrows 
et al. 2000; Bodenheimer, Lin, & Mardhng 2001; Bo- 
denheimer, Laughlin, & Lin 2003; Burrows, Sudarsky, & 
Hubbard 2003; Baraffc et al. 2003; Burrows et al. 2004; 
Chabrier et al. 2004; Laughlin et al. 2005; Baraffe et al. 
2005). 

^ Department of Astronomy and Steward Observa- 
tory, The University of Arizona, Tucson, AZ 85721; 
burrows@zenith.as.arizona.edu, budaj@as.arizona.edu, 
hubeny@aegis. as .arizona.edu 

^ Astronomical Institute, 05960 Tatranska Lomnica, Slovak 
Republic 

Department of Planetary Sciences and the Lunar and Plane- 
tary Laboratory, The University of Arizona, Tucson, AZ 85721; 
hubbard@lpl. arizona. edu 

^ see J. Schneider's Extrasolar Planet Encyclopae- 
dia at |http: / /www. obspm.fr/encycl/encycl. html and the 
Carnegie/California compilation at http;//exoplanets.eu 

^ We will not use the term "Hot Jupiters," since it is a mis- 
nomer from the point of view of spectra and atmospheres. The hall- 
marks of Jupiter's atmosphere, cold ammonia and water clouds and 
methane gas, could not survive in the extreme irradiation regime 
of a close-in EGP that experiences ~10^ times higher fluxes and 
has atmospheric temperatures that are an order of magnitude hot- 
ter. In the atmospheres of the known transiting EGPs, carbon 
is in carbon monoxide, alkali metals, not seen in Jupiter, are the 
predominant absorbers at optical wavelengths, and water is steam. 



Recently, observers and theorists alike have focussed 
on the apparent discrepancy with published theory of the 
transit radii of some EGPs, notably HD209458b, HAT- 
P-lb, and WASP-lb (Knutson et al. 2006; Bakos et al. 
2006b; Charbonneau et al. 2006c), i.e. that these close-in 
EGPs are larger than most theories would predict. Many 
would explain this anomaly by invoking an extra heat 
source for the interior, perhaps caused by orbital tidal 
forcing (Bodenheimer, Laughlin, & Lin 2003), obliquity 
tides when in a Cassini state (Winn & Holman 2005), 
or penetration of gravity waves into the planetary inte- 
rior that then dissipate at depth (Guillot & Showman 
2002; Showman & Guillot 2002). Such a power source 
could indeed be operative, and the powers required are 
not large (|J7|). However, the transit radius of an EGP 
depends upon Mp, the stellar flux at the planet (Fp, 0, 
its atmospheric composition (^, the possible presence 
of an inner core (S}6|), its age, and the atmospheric circu- 
lation that couples the day and the night sides (S}8]). It 
also depends upon the fact that the transit line of sight 
cuts the chord of the planet, and not its radial profile 
O. This effect can add ~3% to ~10% to the measured 
radius (Burrows, Sudarsky, & Hubbard 2003; Burrows et 
al. 2004; Baraffe et al. 2003) and should be included in 
a detailed comparison with observation. 

With so many determinants of a planet's radius, com- 
parison between theory and measurement must be multi- 
parametric. Furthermore, errors in the measured Rp and 
Mp, in the ages, and in the stellar metallicities can be 
large. These introduce significant noise in the interpre- 
tation of any one transiting EGP and are why it is more 
fruitful to look broadly at the entire family. In this way, 
we are able to determine the overall systematics in the 
structures of close-in EGPs and discover trends and char- 
acteristics that would otherwise be obscured if we had 
focussed on one object at a time. As a result, we put 
less weight on our object-by-object "best-fits" than in 
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the patterns that emerge from our study of them collec- 
tively. 

We find that the range of observed radii for the en- 
tire cohort of transiting EGPs is too large to accommo- 
date only one radius anomaly. We show, in fact, that 
some (most) transiting EGPs are smaller than past the- 
ory would have predicted, while we confirm that some are 
larger than past theory would have predicted. We can 
explain both anomalies with 1) enhanced atmospheric 
opacities for the larger EGPs and 2) "ice/rock" cores 
for the smaller EGPs. Such cores are predicted by the 
"core-accretion" model of giant planet formation (Pol- 
lack et al. 1996) and ice/rock cores shrink an EGP of 
a given total mass monotonically with core mass. An 
extreme case is HD149026b (Charbonneau et al. 2006b; 
Fortney et al. 2006). Interestingly, we derive a rough cor- 
relation between the inferred core masses and the parent 
star metallicity. This trend suggests their origin. 

Larger atmospheric abundances, such as those mea- 
sured for Jupiter and Saturn (fJS) Atreya et al. 2003; 
Atreya 2006; Flasar et al. 2005), would lead naturally 
to larger atmospheric opacities that retard the loss of 
heat and entropy from an EGP and delay the shrink- 
age of its radius. However, in this paper we are not 
tying such enhanced atmospheric opacities solely to en- 
hanced atmospheric abundances/metallicities. This is an 
important point. Rather, we are suggesting that the at- 
mospheres of close-in EGPs could also be altered signif- 
icantly by strong optical and UV irradiation. The thick 
hazes, absorbing clouds, and non-equilibrium chemical 
species that could thereby be produced might lead to 
significant increases in the optical thickness of the atmo- 
spheric blanket, leading to a slowdown in the rate of loss 
of core heat. Enhanced atmospheric opacity has an effect 
similar to extra core power. We note that the UV flux at 
the surface of the transiting planets is as much as a factor 
of lO"' times higher than that at the surface of Jupiter. 
Despite the much lower Jovian UV insolation, its atmo- 
sphere contains as-yet-unidentified trace non-equilibrium 
species at the part in '^10^" level that nevertheless result 
in a decrease by almost a factor of two in its blue and 
green geometric albedos. What might be the response 
in the atmosphere of a close-in EGP to the factor of 10* 
increase in UV irradiation? 

Therefore, in this paper we explore the consequences 
for EGP radii of enhanced atmospheric opacities. We 
do this by calculating models using solar, 3 x solar, and 
10 X solar abundance atmospheres, but the latter two 
should be considered ersatze for the effects of enhanced 
opacities of whatever origin. Hence, we decouple the 
effects of increased atmospheric opacity from increased 
envelope heavy-element abundances. If the increase in 
atmospheric opacity were due solely to increased metal- 
licity and our equilibrium chemistry and opacity algo- 
rithms were correct, then the implied increases in the 
heavy-element burden of the envelope, if the heavy frac- 
tion in both atmosphere and envelope were the same, 
could partially or wholly cancel the expansion effect of 
enhanced atmospheric opacity (see ^Hland Fig. [5]). We 
leave open the detailed reasons for the enhanced opaci- 
ties, which could, in addition to super-solar metallicities 
in the atmosphere, be non-equilibrium chemistry, errors 
in the default opacities, and/or thick hazes or clouds. In 
the near future, measurements of both the reflected light 



and thermal emission of close-in EGPs should help to 
constrain both the opacities and compositions of their at- 
mospheres. We note that the detection by Charbonneau 
et al. (2002) of sodium in the atmosphere of HD209458b 
is best fit by the presence of hazes (Fortney et al. 2003) 
or some additional grayish absorber. The default theory 
using clear atmospheres does not explain the factor-of- 
three discrepancy (from merely solar!) in the inferred 
abundance of sodium in HD209458b's atmosphere. 

The upshot of these dual themes concerning atmo- 
spheres and cores is a theory that might explain all the 
transit radii without resorting to an extra power source 
to inflate them. Though an extra power source is still 
possible, we find no simple correlation between the mag- 
nitude of the needed power and any planetary or stellar 
properties. 

In 121 we review the transit data and summarize their 
interesting features, particularly those that demand spe- 
cial explanation. Section [3] demonstrates the general de- 
pendence of transit radii on Mp and stellar flux (Fp). 
The latter varies by more than an order of magnitude 
among the known transiting EGPs. In 21 we discuss the 
"transit-radius" effect that arises from the fact that we 
measure an impact parameter and not a radius. In iJ51 we 
present the results of our calculations without cores for 
solar-opacity and 10 x solar-opacity atmospheres. These 
models are the baseline suite that set the stage for the 
discussions that follow ^ . The higher opacity models can 
flt the large-radius EGPs, modulo remaining uncertain- 
ties in their ages. We note again that we use increased 
atmospheric metallicity as a convenient substitute for en- 
hanced opacity. In ^ we discuss the effects of a cen- 
tral "ice/rock" core and calculate a range of core masses 
needed to achieve better fits for the relevant EGPs. This 
section motivates a possible correlation between the in- 
ferred core masses and the stellar metallicity that might 
inform models of their formation. In fjTl we discuss the 
possible effect on planet structure of an extra heat source 
and determine how much power, object by object, would 
be needed to explain the measured RpS for simple models 
of planet cooling. This section is meant merely to pro- 
vide the reader with a gauge of the range of powers that 
might be required should our default and preferred set 
of models be shown in the future to fail in some crucial 
particular. Curiously, we find that inner cores are still 
suggested by the data even when an extra internal heat 
source is present. In fJHl we discuss a major theoretical 
uncertainty - the advection of heat from the day to the 
night sides due to global circulation. Atmospheric winds 
at altitude and at depth remain wild cards in the general 
theory of EGPs. In SjQl we summarize our results and 
conclusions and reiterate the remaining caveats concern- 
ing the theory of EGP radii. 

2. MEASUREMENTS OF CLOSE-IN GIANT PLANETS 

Table 1 is a compilation of relevant data for the four- 
teen known transiting planets, listed in order of increas- 
ing semi-major axis. These data include semi-major axis 
(a), period (P), Mp, Rp, Fp, and recent observational 
references. We also provide the latest error bars for Mp 

^ We have also calculated 3 x solar-opacity models to better de- 
termine the opacity dependence of transit radii and provide a more 
comprehensive view, but do not provide the corresponding plots. 
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and Rp, though when it seemed prudent we have rounded 
both these and the central estimates. Note that the flux 
at the planet is not monotonia with orbital distance, re- 
flecting the fact that these EGPs orbit a variety of stars 
with luminosities that span an order of magnitude. Ta- 
ble 2 provides these luminosities (L,), along with other 
useful stellar parameters, such as spectral type, stellar 
radius (R*), effective temperature (Tcff ), surface gravity 
(g), metallicity ([Fe/H]), and stellar mass (M). We also 
provide in Table 2 the system distances; some (such as 
those for the OGLE set) are quite approximate. Most 
of the data in Tables 1 and 2 are necessary to construct 
theoretical models and compare them with the measure- 
ments. For instance, to incorporate the effects of stellar 
irradiation one needs models of the stellar spectra and 
luminosities. We employed those of Kurucz (1994), or 
generated our own using the atmosphere code TLUSTY 
(Hubeny & Lanz 1995). In Table 2, we include ages and 
their error bars, both of which should be considered very 
approximate. We list only the central guesses of the stel- 
lar metallicities given in the literature, but ample error 
bars for them should also be assumed. The ages and the 
metallicities are the least well-known quantities in Table 
2 and ambiguities in them translate into uncertainties in 
the interpretation of the theoretical models and transit 
data for any given object. However, as Table 2 suggests, 
the metallicities of these EGP parents vary by a factor 
of ~4. The ages probably range even more broadly. 

Figure [1] depicts Rp versus Mp for all the transiting 
EGPs given in Table 1, along with error bars. Jupiter 
and Saturn arc included for context. This figure encapsu- 
lates the basic measurements to be explained by theory 
and warrants some discussion. The first thing to note 
is that the spread in transit radii is wide, ^40% for the 
bulk and approximately a factor of two when HD149026b 
is included. Some have noted that there is a tendency 
for the larger EGPs to be orbiting the more massive pri- 
mary stars (see Table 2). This is most easily explained 
by the fact that such stars have higher luminosities, and, 
hence, that their EGPs find themselves in more intense 
irradiation regimes (all else being equal) , but planet /star 
distance and planet mass also play central roles. In fact, 
the largest values of Fp are for OGLE-TR-56b, OGLE- 
TR-132b, and WASP-lb, while HD209458b and TrES-2 
are in the middle of the pack (see Table 1). As the upper 
envelope of the data in Fig. [T] suggests, there is a slight 
tendency for the lower mass EGPs to have higher radii. 
This effect is a straightforward consequence of basic the- 
ory and is at least as important 

There are other apparent curiosities. Using Fig. [T]and 
Table 1, we can compare subsets of EGPs with roughly 
the same Mp. One such triplet, in order of decreasing 
radius, is WASP-lb, XO-lb, and WASP-2b. We might 
expect that, given this radius hierarchy, Fp would mono- 
tonically decrease from WASP-lb to WASP-2b. How- 
ever, Fp for XO-lb is lower than that for WASP-2b. 
HAT-P-lb, OGLE-TR-lOb, and OGLE-TR-lllb consti- 
tute a similar triplet, but Fp for OGLE-TR-lllb is the 
largest of the three, breaking what should otherwise be 
a monotonia trend. Moreover, the radii and masses of 
HD189733b and OGLE-TR-132b are roughly the same, 
yet their FpS are almost an order of magnitude differ- 
ent. The most extreme case is HD149026b, which has the 
fourth highest Fp, but the smallest radius. Our overall 



thesis is that these features can be explained, to within 
the error bars, not one-dimensionally, but only after the 
various effects of Mp, Fp, core mass, atmospheric opacity, 
and age are simultaneously addressed. 

Figure [5] depicts the dependence of the measured Rp 
on the estimates of the stellar metallicity ([Fe/H], Table 
2). Error bars in both quantities, in particular [Fe/H], 
will smear this plot, but the basic relationships, if there 
are any, should emerge as plotted. We see that at all 
metallicities, there is a wide range of measured radii, 
and no clear and simple correlation with either Mp or 
Fp (Table 1). Curiously, there seem to be two branches 
(upper and lower), but this may be an artefact of small- 
number statistics. In any case. Figs. [1] and [Hand Tables 
1 and 2 summarize the salient information concerning 
the known transiting EGPs to be explained by theory. 

3. DEPENDENCE ON STELLAR FLUX AND Mp 

To demonstrate the general dependence of Rp upon or- 
bital distance and Mp, we have generated Fig. [3l In it, we 
depict evolutionary trajectories for a Saturn-mass planet 
(0.3 Mj ^, solid) and a Jupiter-mass planet (dashed) at 
distances from a G2V main sequence star of 0.02, 0.03, 
0.04, 0.05, and 0.06 AU. These models are not per se 
our preferred models for any of the known transiting 
EGPs, assume solar-metallicity atmospheric abundances 
(Asplund, Grevesse, & Sauval 2006) and opacities, do 
not include inner cores, but, as do all the models we 
present in this paper, employ the well-developed bound- 
ary condition formalism of Burrows, Sudarsky, & Hub- 
bard (2003) and Burrows et al. (2004). For these, and aU 
evolutionary calculations in this paper, we pre-calculate 
grids of self-consistent irradiation boundary conditions at 
130 points that span the internal flux and surface grav- 
ity space (Burrows, Sudarsky, & Hubbard 2003) likely to 
be traversed during the evolution of each single primary 
star/semi-major-axis combination. During each evolu- 
tionary calculation, we interpolate in this grid of bound- 
ary conditions. Appropriately different stellar spectra 
(see Table 2; Kurucz 1994; Hubeny & Lanz 1995) for 
each system are employed and we set up these grids for 
each of the 14 known transiting EGPs and three sets of 
opacities (|J5]). Hence, for this study we have calculated 
14x130x3 — 5460 detailed spectral/ atmosphere models. 

We see immediately that the radius of a low-mass EGP 
is more sensitive to distance, with that of a Saturn-mass 
EGP varying by ~0.2 Rj » from 0.02 AU to 0.06 AU 
and that of a more-massive Jupiter-mass EGP varying 
by '-^0.1 Rj over the same orbital distance range. More- 
over, younger EGPs have larger radii than older repre- 
sentatives, but after '^1.0 Gyr all evolutionary trajecto- 
ries start to flatten. This fact emphasizes the potential 
role of youth in providing large radii, and the ambigu- 
ities that arise in the interpretation of transiting EGPs 
with poorly-known ages. This is particularly relevant for 
OGLE-TR-lllb, HD189733b, TrES-2, WASP-lb, and 
WASP-2b, whose ages are either unknown or very poorly 
known. Figure [3] also shows that the timescale for radius 
decay is longer for lower-mass EGPs. 

Figure [4] continues our demonstration of the effects of 
irradiation and planet mass on Rp by depicting its direct 

1 Mj (Jupiter's mass) = 1.89914x10^" g 
* 1 Rj (the radius of Jupiter) = 7.15x10^ cm 
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dependence on the stellar flux (Fp) at the substellar point 
for the same class of theoretical models. Roughly one 
order of magnitude in Fp is depicted. Masses of 0.3, 
0.5, 0.65, 1.0, and 1.25 Mj are shown for an age of 2.5 
Gyr. This age is roughly the mean age of stars in the 
solar neighborhood. Again, we see that, all else being 
equal, smaller-mass EGPs have larger radii and depend 
more steeply upon Fp. For a 1.25-Mj EOF, Rp varies for 
the depicted range of FpS by '^0.08 Rj , while for a 0.3- 
Mj EGP it varies by as much as ^0.24 Rj . This behavior 
is consonant with our statement in ij2] that the upper 
envelope of the data depicted in Fig. [T] has a negative 
slope. Note that the spread in Rp with mass at high Fp 
is significantly larger than at low Fp. This is connected 
with the convergence of the radii of cold EGPs due to 
the n = 1 polytropic character of the H2/He equation of 
state (Burrows et al. 2001). 

4. TRANSIT RADIUS EFFECT 

Measuring a transit provides the impact parameter of 
the planet, not its photospheric radius. This means that 
the planetary limb, through which the light from the star 
that defines the depth of the transit emerges, is at a 
slightly larger distance from the projected planet cen- 
ter than the canonical r = 2/3 planetary radius. For 
large FpS, high atmospheric metallicities, and small Mp, 
this difference can be ^5%. Hence, the effect should be 
included in any comparison with data and failure to in- 
clude it will exaggerate the apparent discrepancy with 
the previous theory of the radii of the biggest transiting 
EGPs. For the 0.64-Mj EGP HD209458b, the effect can 
be larger than 0.05 Rj (Burrows, Sudarsky, & Hubbard 
2003; Baraffe et al. 2003). 

The wavelength-dependent optical depth, Tchord, along 
a chord followed by the stellar beam through the planet's 
upper atmosphere, is approximately: 



Tchord ~ KPphH^ ^ " ' ^^"^ 

where k is the wavelength-dependent opacity, pph is the 
mass density at the photosphere, ARch is the excess ra- 
dius over and above the Tph — | radius (the radius of 
the traditional photosphere), and H is the atmospheric 
density scale height. The latter is given approximately 
by kT / pgrup, where fi is the mean molecular weight, g is 
the surface gravity, T is some representative atmospheric 
temperature, and nip is the proton mass. By definition, 
and ctssumiiig a.11 GxponenticLl atmosphere, Tpfi — ^Pph 

H 

= |. For Tchord to equal |, this yields 



AR,h^Hln^—^^5H. (2) 

ARch should be included in the theoretical radius that is 
compared with the measured transit radius. In this pa- 
per, we include it implicitly by first calculating the radius 
of the convective-radiative boundary and then adding 
to it the additional distance to the Tchord — 2/3 level 
in the corresponding detailed atmosphere model. We 
refer to this additional distance as AR (no subscript), 
which contains AR^h- Figure [H] depicts AR versus planet 
mass (Table 1) for solar (black) and lOxsolar (red) at- 
mospheric opacities and representative coreless models 



of twelve of the measured transiting EGPs. The dis- 
tance, AR, from the radiative-convective boundary to 
the Tchord = 2/3 level is between ~0.04 Rj and ~0.15 
Rj for H2 /He-dominated atmospheres, depending mostly 
on the planet's mass {Mp), the stellar flux at the planet 
(Fp), and (weakly) its age. As Fig. [5] indicates, AR is 
smaller for higher-mass EGPs, larger for planets experi- 
encing higher FpS (see Table 1), and larger for higher at- 
mospheric opacities. Concerning the latter, the increase 
in AR in going from solar to lOxsolar ranges from ~0.01 
to ^0.04 Rj . The contribution of ARch to AR varies 
from ^^10% to ^50%. Note that the numbers depicted 
in Fig. [5] assume that the mean molecular weight (fi) is 
not altered at high opacity. Even if high opacity meant 
high metallicity, the p effect at lOxsolar would amount 
to a diminution of the scale height and the transit-radius 
effect itself by no more than ~20% of the enhancement, 
and would not compensate for the corresponding increase 
in AR due to the opacity effect. 

5. MODELS WITH SOLAR- AND ENHANCED-OPACITY 
ATMOSPHERES AND NO CORES 

In situ and remote-sensing measurements of the at- 
mospheric compositions of the giant planets Jupiter and 
Saturn reveal that most of the dominant elements, such 
as carbon, nitrogen, and sulfur, exist there in super-solar 
abundances (Atreya et al. 2003). Atreya (2006) esti- 
mates that [N/H] and [C/H] in Jupiter's atmosphere are 
4-5 times solar and that [C/H] in Saturn's atmosphere 
is 9-10 times solar. Furthermore, Flasar et al. (2005) 
estimate that carbon in Saturn's atmosphere is ~7 times 
solar. Given the ambiguities in the interpretation of the 
Galileo probe results, [0/H] is problematic, but it too is 
widely considered to be super-solar. Since the metallici- 
ties of EGP host stars are preferentially in excess of the 
Sun's (Fischer & Valenti 2005), the idea that the atmo- 
spheres of orbiting EGPs are heavy-element-rich is more 
than just an intriguing possibility. In addition, the ex- 
cesses seen in Jupiter and Saturn are in keeping with the 
core-accretion model of giant planet formation (Pollack 
et al. 1996), and are some of the reasons it is preferred. 

It was these super-solar heavy-element abundances in 
the Jovian planets that first motiviated us to explore 
the effects on EGP radii of enhanced atmospheric opac- 
ities. As we suggest in Sjl] even for solar abundances, 
strong irradiation may significantly alter the chemistry 
and opacities of the atmospheres of close-in EGPs. Here- 
after, we use super-solar metallicity as a substitute for 
enhanced opacity by whatever means and for whatever 
elemental abundance pattern and metallicity. We ex- 
plore the consequences for the radii of irradiated EGPs 
of such opacities and compare with the corresponding 
results for default solar-metallicity atmospheres. In the 
models that follow, 3 x solar and lOxsolar are to mean 
"with heavy-element opacities that are 3 and 10 times 
what they would be at a given temperature and pres- 
sure for the canonical, unaltered solar-metallicity atmo- 
sphere." However, note that the envelopes of the models 
presented here are assumed to be pure H/He mixtures 
and that the effect on the planet's radius of envelope 
metals is, for our purposes, "absorbed" into an effect 
due to the core alone. Hence, our cores "stand in" for 
the core/envelope vis a vis their summed effect on the 
planet radius (Sj6]). 
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Higher atmospheric opacities retain the core's heat and 
entropy, and this maintains the EGP's radius at higher 
values for longer times. This consequence of higher at- 
mospheric gas-phase opacities (which could be abetted 
by upper atmosphere clouds; cf. Fortney et al. 2003) 
is similar in effect to that of an extra core power source 
Oi but we believe that this explanation of large EGP 
radii may be more natural. 

As stated in for all our calculations, we em- 
ploy the evolutionary, spectral, atmospheric, and opacity 
techniques described in Burrows, Sudarsky, & Hubbard 
(2003) and Hubeny, Burrows, & Sudarsky (2003), and 
discussed in Burrows et al. (2001)^. We set the redistri- 
bution factor (Burrows et al. 2004), /, equal to 1/4, and, 
therefore, assume complete heat redistribution at depth 
(see ^ . Figures [5] and [7] portray theoretical evolution- 
ary trajectories of Rp versus age for coreless models of 
all fourteen of the known transiting EGPs. Model atmo- 
spheres for both solar opacity (top) and 10 x solar opacity 
(bottom) are shown and models with (right) and with- 
out (left) the AR term are included for comparison. The 
measured transit radii and ages are superposed, along 
with error bars (Tables 1 and 2). For each EGP, the 
color used for both model and data is the same. Since 
the ages for WASP- lb and WASP-2b are not in the lit- 
erature, we arbitrarily set them equal to 2 ± 1 Gyrs. 

Figure[S]contains eight of the smallest measured EGPs, 
and Fig. [7]contains the other six (and, hence, the largest) 
EGPs. The cut between the two sets is of no fundamental 
significance. As Fig. [^indicates, if we use solar opacities, 
ignore AR, and leave out a core, the left-hand-side mod- 
els would fit the corresponding data rather well, except 
for HD149026b. All the coreless models of HD149026b 
are discrepant by wide margins (by as much as a factor 
of two) and a core of substantial mass seems the only 
option (Fortney et al. 2006; ©. In fact, HD149026b is 
more like a super-Neptune than an EGP. 

The wide range of possible ages for some of the EGPs 
depicted in Fig. [6l in particular for OGLE-TR-lllb and 
0GLE-TR-113b, makes interpretation a bit uncertain, 
particularly in the lower age range. However, for longer 
ages the models are substantially age- independent. In 
addition, one can't arbitrarily ignore the AR term and 
as the top right-hand panel of Fig. [6] indicates, the 
solar /coreless "fits" then evaporate when including it. 
Even if the errors in Rp are obliging, and data and model 
for a few of the eight EGPs are reconciled, one is unlikely 
to be able to do this for all of them. The upshot is that 
even at solar opacities coreless models for these smaller 
transiting EGPs are disfavored. Models portrayed in the 
10 X solar panels at the bottom of Fig. [6] are even more 
disfavored. The actual opacities of the atmospheres of 
these EGPs don't have to be as high as for the 10 x solar 
models for these plots to be indicative of a severe prob- 
lem. This is the first major radius problem: many of the 
known transiting EGPs are too small, not too large. 

Figure [7] depicts the six largest transiting EGPs in the 
same format as Fig. |6l The gap with theory for so- 
lar opacities, no cores, and no AR term is wide for all, 
except for TrES-2, if its age is quite low. As the up- 
per right-hand panel indicates, including the AR effect 

^ We assume that the stellar luminosity does not evolve with 
time. 



helps, but not enough. However, at '--^lOxsolar, the mod- 
els for all these larger EGPs start to fit rather well, the 
degree of fit depending centrally upon the age and radius 
error bars. In fact, for OGLE-TR-lOb and OGLE-TR- 
56b, their 10 x solar-opacity radii are on average too large. 
This is true for OGLE-TR-56b, despite its large Fp (Ta- 
ble 1). The measured radius of HD209458b is stiU a bit 
larger than the theory, but it is within 1.5-cr for its cen- 
tral age estimate, and better than this for younger ages. 
HAT-P-lb fits well, TrES-2 fits well for a wide range of 
ages. WASP-lb can fit, in particular if it is not very 
old (recall that its age is unknown). The largest opac- 
ity effects, those associated with an increase in radius of 
~0. 05-0.1 Rj in going from solar to lOxsolar, obtain for 
the least massive EGPs (lowest Mps) with the highest 
irradiation fluxes, FpS (see Table 1). For HD149026b, 
for which Mp =0.36 Mj and Fp is the fourth highest, 
the magnitude of the atmospheric opacity enhancement 
effect is ~0.2 Rj . 

Therefore, we conclude that higher-opacity atmo- 
spheres and the inclusion of the AR term can explain 
the largest of the measured radii. Moveover, the range 
of radii among the fourteen known transiting EGPs is too 
wide to be explained by one factor alone. Importantly, 
there is a small-radius problem as well, one that can not 
be solved by an extra heat source. We next show in fj6] 
that ice- rock cores for almost all the known EGPs, with 
smaller cores for the largest EGPs, are indicated. 

6. EFFECT OF A CENTRAL CORE 

In the core-accretion model of giant planet formation 
(Pollack et al. 1996), a mass of ice and rock accummu- 
lates until it achieves a critical mass. This critical mass 
then nucleates rapid gas accretion and the giant planet 
grows to its final mass at the expense of the surrounding 
protostellar nebula. Such a two-step process is suggested 
because nebular temperatures are estimated to be too 
high for the inferred disk areal mass densities to allow 
direct gravitational instability by the Toomrc condition 
(Boss 1997, 2001), akin to the Jeans criterion for star 
formation. Importantly, there is direct evidence for the 
presence of a ~15 Earth-mass^° core in Saturn and some 
evidence for a similar core in Jupiter (Guillot & Saumon 
2004). The ice giants Neptune (17.1 Earth masses) and 
Uranus (14.5 Earth masses) are thought to be such nu- 
clei that may have been starved of gas at birth by the 
low-density neighborhood in which they were born. 

In all cases, for a given total planet mass, the presence 
of a core shrinks the total radius of an EGP. We numer- 
ically incorporate such cores into our models by placing 
a compressible ball of olivine in the center of the model 
planet. For each model, the core mass (Mc, in Earth 
masses, as per convention) is set and pressure continuity 
between the solid core and the gaseous envelope is en- 
sured throughout the evolution. The ANEOS equation of 
state (Thompson & Lauson 1972) is used for olivine and 
the Saumon, Chabrier, & Van Horn (1995, SCVH) equa- 
tion of state is used for the H2/He envelope. In these cal- 
culations, we assume that the specific heat capacity per 
mole of the solid cores is the same as derived using the 
SCVH equation of state. What the actual specific heats 
and entropies of the core and heavy-element component 

One Earth mass is equal to 5.98x10^^ g. 
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of the envelope are is an important open issue. If the core 
has a high thermal inertia, this can delay the cooling of 
the planet and the shrinkage of its radius. Conversely, 
if the heat capacity of the core is smaller than that of 
H/He mixtures, large core models will cool down slightly 
more quickly than our corresponding models, resulting in 
slightly smaller planet radii. The zero-pressure density 
of olivine is ~3.2 g cm~^, significantly higher than the 
average density of EGPs (Charbonneau et al. 2006a). 
This is the point. If we replace the olivine with ices or 
ice/rock mixtures the results vary slightly, but not qual- 
itatively. Reliable equations of state for heavy-element- 
rich gaseous envelopes that could constitute most of the 
planet's mass are still not available, so we assume that 
these envelopes are dominated by H2/He mixtures. We 
have set the helium mass fraction equal to 0.25. Some 
think that whether the heavy elements are in the core 
or the envelope, their effect on Rp is the same. This 
has not been shown, but one can consider the inner core 
masses with which we deal as substitutes for the total 
heavy-element burden in the planet. It is the systemat- 
ica in the group of known transiting planets for which we 
are looking and the favored parameters of each EGP are 
bound to improve significantly with time. 

Figure [8] plots theoretical total radii as a function of 
core mass, Mc, for the estimated ages of OGLE-TR- 
10b, OGLE-TR-56b, HD189733b, and XO-lb (Table 2). 
These are merely representative. The lines in Fig. [5] are 
for solar, 3 x solar, and 10 x solar models. The measured 
radii of these transiting EGPs are given as dots and the 1- 
cr radius error bars are indicated with vertical lines. The 
dots are placed arbitrarily along the horizontal direction 
at core masses equal to the mass fraction represented by 
3 X solar metallicity times the total EGP mass and the 
rightmost extent of the horizontal "error bars" is placed 
at the corresponding 3 x sie^/ar metallicity masses. If the 
central value of the estimated stellar metallicity is below 
solar (e.g., HD189733b), the horizontal line is truncated 
at the dot. As Fig. [^indicates, a core mass of ^20 Earth 
masses can shrink an EGP by ~0.05— 0.1 Rj . EGPs 
with smaller total masses (such as for OGLE-TR-lOb 
and HD149026b) manifest a steeper drop in radius with 
increasing M^. Also, while only small cores are indicated 
for XO-lb and HD189733b, for OGLE-TR-lOb the core 
mass derived for an atmospheric metallicity of ~3x solar 
(if metallicity and opacity were tied) is also the preferred 
core mass, i.e., the center of the yellow cross intersects 
the dashed line. This is not often the case in our current 
model set, since we have decoupled envelope metallicity 
from atmospheric opacity. However, note that, as Fig. 
[8] shows, since the intercept of the solar-metallicity lines 
with the y(radius)-axis is often (though not always) be- 
low the radius positions for, e.g., the 3 x solar metallicity 
dots for the dashed 3 x solar lines, the increase in the ra- 
dius due to increasing the metallicity in the atmosphere 
is often slightly larger than the decrease in the radius due 
to a possible corresponding increase in the metallicity in 
the envelope. Moreover, adding ices to constitute a true 
"ice/rock" core, substituting for the pure olivine core we 
now assume, should slightly favor larger radii, but we 
leave this to future studies. So, higher metallicities over- 
all can still be an important part of the solution to the 
large-radius problem. Nevertheless, more work on the 
envelope equation of state for arbitrary heavy-element 



fractions is still clearly needed. 

Table 3 lists in bold the approximate core masses that 
provide model fits with solar, 3 x solar, and 10 x solar at- 
mospheres for each of the fourteen EGPs. We have 
rounded the best-fit core masses to the nearest conve- 
nient number. In parentheses in each column, to the left 
and the right of the bolded values, are the best-fit core 
masses for radii that are il-a from the central radius es- 
timate (Table 1). Hence, the right value is for the larger 
(-|-l-cr) radius and the left value is for the smaller (-1-cr) 
radius. If there is no good fit, or Mc would have had to 
be negative, we print (••■). "0" means close to zero, but 
could be a bit larger. Since we have no age estimates 
for WASP-lb and WASP-2b, we leave the corresponding 
rows for them empty in anticipation of future data. For 
HD149026b, we provide only central estimates and for 
HD209458b (for which we provide no estimates) the best 
fits require that the actual transit radius be beyond its 1- 
tr radius error bars. All these core masses are derived at 
the central age estimates given in Table 2. Since in most 
cases these ages are quite uncertain, the actual ages could 
yield very different core mass estimates. For instance, if 
an EGP's age is significantly younger, the predicted ra- 
dius without a core would be higher (see Figs. [S] and 
[7]). In that case, compensating for the resulting larger 
radius deficit would require a larger core mass, all else 
being equal. 

We see in Table 3 that larger core masses are required 
in models with higher atmospheric opacities, with a 
swing of ^20-30 Earth masses from solar to 10 x solar. We 
also see that the range of theoretical values for Mc is very 
large, from zero to ^^100 Earth masses. Furthermore, 
Table 3 suggests that the canonical "15-Earth masses" 
that works for our solar system giants might be disfa- 
vored as the generic giant planet core mass. Moreover, 
we note that the high-Mp OGLE-TR-132b and the low- 
Mp HD149026b both require very large cores, though su- 
perficially the radius of OGLE-TR-132b might not have 
seemed anomalous. For HD149026b, with a small Mp, 
a large Fp, and a small Rp, the conclusion that a large 
core is required is unexceptional. But for the more mas- 
sive OGLE-TR-132b, with the highest Fp of the family, 
it is intriguing that a very large Mc of comparable mag- 
nitude may be required. We draw a similar conclusion 
for 0GLE-TR-113b, which is the most massive of the set 
and has a modest Fp, but may require a core mass of 
60-80 Earth masses. 

What patterns emerge from this theoretical study and 
Table 3? Figure [9] plots the parent stellar metallicity 
versus the theoretical core masses given in Table 3 for 
twelve of the known transiting EGPs. The different dots 
for each planet are for the three different atmospheric 
opacities. In this plot and in Table 3, the dependence 
of Mc on atmospheric opacity for each of the EGPs is 
seen to be less important than the wide spread in Mc 
from object to object. The most important feature to 
emerge from Fig. [9] is that Mc seems to increase with 
[Fe/H]. Based upon their preliminary analysis, Guillot et 
al. (2006) suggest a similar correlation. Those stars with 
the lowest [Fe/H], such as HD189733, XO-1, HD209458, 
and the parent of TrES-2, all seem to be orbited by EGPs 
that require small cores. Those stars with the largest 
values of [Fe/H], such as OGLE-TR-132 and HD149026, 
seem to house EGPs that require the largest cores. A 
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"straight" line can be drawn through the points, suggest- 
ing a correlation between inferred core mass and stellar 
metallicity. Moreover, around solar values of the stellar 
metallicity, the suggested core masses are in the solar- 
system regime, paralleling Jupiter and Saturn. Finally, 
on Fig. [9] at low stellar metallicity no large cores are de- 
rived, and at high stellar metallicity no small cores are 
derived. To be sure, there are deviations from this sim- 
ple picture, such as OGLE-TR-56b and OGLE-TR-lOb, 
but these points are derived using central values of the 
poorly known ages and stellar metallicities. If the metal- 
hcities and/or ages of OGLE-TR-56 and OGLE-TR-10 
are slightly lower, the corresponding points will move up 
and to the left, into the trend line. Similarly, if we derive 
their core masses using the upper 1-cr radii, the corre- 
sponding dots will shift upward on Fig. [S) However, it is 
also not altogether unreasonable to expect some scatter 
in giant planet formation and in Mc. 

On Fig. [9l for each single object one point separately 
is not very suggestive, but plotted together they collec- 
tively indicate a correlation that hints at their origin. 
At the very least, super-solar and super-stellar heavy- 
element abundances in the interiors of these planets, if 
not the presence of cores per se, are strongly suggested. 
Hence, we offer Fig. |9] as tantalizing evidence for the 
presence of dense cores and/or heavy-element-rich en- 
velopes in EGPs and, therefore, for the core-accretion 
model of giant planet formation. 

7. EFFECTS OF EXTRA HEAT SOURCE IN INTERIOR 

Many workers have sought to explain the large radii 
of transiting planets such as HAT-Plb, WASP-lb, and 
HD209458b by invoking an extra power source in the 
planet's interior (Bodcnheimer, Laughlin, & Lin 2003; 
Guillot & Showman 2002; Winn & Holman 2005; 
Chabrier et al. 2004; Charbonneau et al. 2006c). Such 
a power source would maintain the entropy in the core, 
and hence its radius, by compensating in part for radia- 
tive cooling at its periphery from all quadrants. As our 
discussions in [J5] and [J6] indicate, we don't prefer this 
solution, and in fact conclude that there are two radius 
problems, only one of which could be resolved with an ex- 
tra core heat source. Nevertheless, it is useful to estimate 
the magnitude of the power required for each transiting 
EGP to affect its measured radius. The goal is to de- 
termine whether the requisite power could be correlated 
with some other system parameter, such as intercepted 
stellar power, Lp. In Table 4, we provide such estimates 
for the transiting EGPs for two cooling models. The first 
(labeled "Power (Iso)") ignores stellar irradiation com- 
pletely and assumes the object can otherwise be consid- 
ered isolated (see also Chabrier et al. 2004). The central 
value of the measured radius (Table 1) is assumed to be 
the target of the fit and AR is not added. Solely for the 
purposes of illustration, the atmospheres have solar opac- 
ities. We see in Table 4 that between 0.45% and 0.005% 
of each EGP's Lp would be called for. The characteristic 
variation is a factor of ten. This needed variation from 
object to object makes unclear the origin of such a power 
source. 

The second model ("Power (Solar)") also assumes that 
the atmospheres have solar composition and drops the 
AR, but includes the effect of stellar irradiation with 
our default redistribution parameter {^). These models 



are the solar-atmosphere/no- AR models described in fj5l 
but with an extra power source. In this case, the range of 
fractions of Lp is more narrow, between 0.01% and 0.05%, 
and a factor of ten smaller than for the "Power(Iso)" 
model set, reflecting the effect of irradiation. Note that 
for more than half the models in this model set an extra 
heat source would make the radius fit worse, not bet- 
ter. Other atmospheric opacities/metallicities could have 
been used in this illustrative study, but the qualitative 
results would have been similar. To further demonstrate 
the dependence on core power of the evolution of Rp, 
Figure \W\ depicts such trajectories for two representive 
EGPs, HD209458b and HAT-Plb, for both "Power (Iso)" 
and "Power (Solar)" assumptions and for a variety of core 
powers. 

While it is noteworthy that the fraction of Lp needed 
to modify Rp in a measureable way is quite small, no 
natural mechanism and no systematic reason for signifi- 
cant variation from object to object suggest themselves. 
Nevertheless, the possibility of an internal power source 
can not yet be eliminated out of hand. Indeed, such ex- 
tra heating may emerge as another degree of freedom in 
the fits. However, at present we find that Occam's Razor 
and the arguments in fJS] and [JHl obviate the necessity for 
a central role for such an ad hoc core power of undeter- 
mined provenance. 

We end this section with a curious observation. On Fig. 
[HI we have placed gold points to indicate the approximate 
core masses necessary to fit models having an extra inter- 
nal power source whose magnitude is an arbitrary, fixed, 
percentage (0.3%) of Lp (the same percentage for all the 
different Lps). These models have solar- metallicity at- 
mospheres, but no irradiation or AR effects. Even for 
these models, we see the same general trend of inferred 
core mass with stellar metallicity that was identified in 

m 

8. AMBIGUITY IN COOLING FROM THE DAY/NIGHT 
SIDES 

The day/night difference in the cooling rates of 
strongly irradiated planets remains the most uncertain 
aspect of all published theories. If there is no modifica- 
tion of the heat flux at the radiative/convective bound- 
ary on the night side due to heat redistribution at depth 
from the day side, and the planet cools on the night side 
as if isolated, then the nightside losses will overwhelm 
the much smaller dayside losses and an extra heat source 
([J7]) may well be required to explain those EGPs with 
the largest transit radii. 

In our default cooling model, we set the redistribution 
parameter, /, defined in Burrows, Sudarsky, & Hubbard 
(2003) and Burrows et al. (2004), and used by other 
groups (e.g., Fortney et al. 2006; Chabrier et al. 2004), 
equal to 1/4. This value signifies complete heat redistri- 
bution at depth and longitude-independent interior core 
fluxes outward. The factor, /, influences the day/night 
temperature (T)) /pressure (P) profile contrasts only at 
high pressures near the radiative/convective boundary 
(at Rosseland rs of ~10^; see Burrows et al. 2004). At 
altitude, the day/night contrast in an EGP's spectrum, 
formed at lower Rosseland ts of 0.1 to a few, can be large, 
as suggested by the recent Ups And b light-curve data 
(Harrington et al. 2006). However, at the same time, 
zonal winds at high optical depths can still efficiently 
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redistribute heat and entropy. It is the T/P profile at 
depth that regulates core cooling. Such efficient deep 
heat transport is suggested in the work of Showman & 
Guillot (2002) and Guillot & Showman (2002), but is by 
no means proven. Nevertheless, we make this assumption 
in order to discover and explain the systematic features 
across the family of known transiting EGPs. 

As an aside, we note that many people think that rota- 
tion is an efficient means to transport heat globally, and 
use Jupiter and Saturn as examples. There is almost 
no latitude or longitude dependence of the mid- and far- 
infrared emissions of either Jupiter or Saturn. Despite 
the secant effect of the incident stellar flux, their emis- 
sion temperatures at these wavelengths are almost com- 
pletely uniform. However, it is not rotation that smooths 
out these emissions, but the direct heating of the convec- 
tive regions of these planets by solar infrared (IngersoU 
1976; Hubbard 1977; and IngersoU and Porco 1978). The 
radiative/convective boundary is at low optical depths in 
these solar-system giants. As a result, the solar heat di- 
rectly absorbed in the convective zone is efficiently re- 
distributed throughout the planet's interior, setting a 
uniform inner boundary for internal heat flux outward. 
For closer-in EGPs, the radiative/convective boundary is 
at greater depths and this mechanism does not operate. 
The upshot is that for more strongly-irradiated EGPs, 
the mechanisms for longitudinal heat transport are more 
subtle, and problematic. A number of groups are at- 
tempting to address this issue with multi-dimensional, 
though approximate, numerical tools (Menou et al. 2003; 
Cho et al. 2003; Burkert et al. 2005; Cooper & Showman 
2005), but these efforts are only in their early stages. 

9. DISCUSSION AND CONCLUSIONS 

In this paper, we have calculated the theoretical evolu- 
tion of the radii of all fourteen of the known transiting gi- 
ant planets for a variety of assumptions concerning their 
atmospheric opacities, inner core masses, and possible 
internal power sources. We have incorporated the effects 
of stellar irradiation and have customized such effects for 
each EGP and star. Using measurements of their ages, 
masses, and transit radii, we have sought to reconcile 
these transit radii with theory. While it can be difficult 
to fit each EGP definitively, looking at them collectively 
can reveal important underlying features of the family 
as a whole. In doing so, we find that there are two, not 
one, radius anomalies. Not only are the radii of a sub- 
set of the known transiting EGPs larger than expected 
from previous theory, but many of the other objects are 
smaller than expected. Unless all the atmospheres have 
only the default ^solar-metallicity opacities, the AR ef- 
fect can be ignored, and an internal power source whose 
magnitude is not correlated in any obvious way with sys- 
tem parameters is operative, we conclude that the spread 
of measured radii is too large not to admit of a dual prob- 
lem. 

We suggest that the larger EGPs can be explained 
by invoking enhanced opacity atmospheres, that might 
be due only in part to enhanced metallicity, that natu- 
rally retain internal heat, and, hence, maintain their radii 
larger, longer. This can be done without an extra inter- 
nal power source, though such a source can not yet be 
eliminated either as an important or a sub-dominant as- 
pect of the theory for some irradiated EGPs. We offer en- 



hanced atmospheric opacities as a more straightforward 
explanation for the large-radius EGPs. Such an expla- 
nation, however, may require non-equilibrium chemistry 
and/or haze formation in the severe irradiation regimes 
in which transiting EGPs find themselves and we have 
not provided in this paper a detailed chemical rationale 
for such altered atmospheres. 

Furthermore, we suggest that the other anomaly, that 
of the small radii we find for the majority of the known 
transiting EGPs, can be explained simply by the pres- 
ence of dense cores and/or metal-rich envelopes in most, 
or all, of these fourteen objects. For no EGP orbiting a 
lower-metallicity star do we infer a large inner core. Con- 
versely, for no EGP orbiting the highest-metallicity stars 
do we infer a small inner core. Moreover, the core masses 
we find for EGPs transiting near-solar metallicity stars 
are close to those estimated for Jupiter and Saturn. Im- 
portantly, we derive a roughly montonically-increasing 
relationship between the stellar metallicity and the es- 
timated core mass. High stellar metallicity has been 
shown to correlate with the probability of the presence 
of an EGP in the radial-velocity data (Fischer & Valcnti 
2005). In this paper, we find that high stellar metallic- 
ity may also imply large inner cores and/or metal-rich 
envelopes. These twin correlations may speak to the 
mechanism of EGP formation and are in keeping with 
the core-accretion model of their origin. 

There are a number of caveats to our conclusions. First 
is the uncertainty concerning the nightside cooling. If 
there is no means by which cooling of the interior can be 
stanched by heat redistribution at depth from the day- 
side (Burrows et al. 2004), then an extra power source 
might be required for the larger radii. Second is the 
wild card of rotation. Since close-in EGPs are no doubt 
in synchronous rotation at periods larger than those of 
Jupiter and Saturn, the effects of rotation will result in 
no more than a few percent expansion, but have not yet 
been included in our analysis. Furthermore, centrifugal 
expansion is most manifest in the transit plane. Third is 
the possibility of delayed migration of some of the plan- 
ets. If migration were to take many tens of millions of 
years (Murray et al. 1998), then the planet might have 
had time to cool and shrink as if in isolation, without the 
benefit of the effects of irradiation. Subsequent irradia- 
tion when in extremis could not reinflate the core (Bur- 
rows et al. 2000). Fourth is the fact that we have merely 
motivated altered chemistry in the atmospheres of these 
severely irradiated EGPs, and have not demonstrated 
the required chemistry, nor the opacity-enhancing effects. 
High-metallicity atmospheres in themselves would be ad- 
equate, but if these were accompanied by envelopes with 
similar metallicities, the radius-increasing effect can be 
partially of wholly cancelled. As Fig. [8] demonstrates, in 
many, though not all, of the cases the enhanced opac- 
ity effect of super-solar metallicity in the atmosphere 
can still trump the shrinkage effect of the same metal- 
licity in the envelope. Super-solar metallicity in the at- 
mosphere, expected generically for EGPs, can still be 
part of the solution to the large-radius problem. How- 
ever, in this study we have decoupled the two and fu- 
ture detailed work on UV-driven chemistry, the opacities 
of strongly irradiated and synchronously-rotating atmo- 
spheres, and the equation of state for general mixtures 
is clearly needed. Fifth is the possibility that the heavy 
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elements and the dominant absorbing compounds of the 
atmosphere might settle gravitationally, thereby deplet- 
ing it of its high-opacity components. Without these 
species, the high-opacity effect that we suggest may be 
instrumental in explaining the largest EGP radii would 
be compromised. However, mixing due to the vigor- 
ous shear motions caused by the zonal winds anticipated 
throughout these regions may in fact be adequate to en- 
sure an unstratified atmosphere. Nevertheless, relevant 
calculations to estimate such mixing are warranted. Fi- 
nally and sixth are the remaining ambiguities in system 
age, EGP radius, and stellar metallicity. The inferred 
core masses, or range of core masses, and the fits to the 
larger-radius EGPs depend upon those parameters. Our 



results could be more robust or less robust, depending 
upon the eventual values of these quantities. 
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Transiting Planet Data'^ 
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TABLE 2 

Data on Parent Stars-'^ 
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TABLE 3 

Approximate Inferred Core Mass Ranges given Central Age Estimates^ 
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TABLE 4 

Internal Power That Would Be Necessary^ 
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Fig. 1. — Transit radii (Rp, in Rj ) of all of the irradiated EGPs listed in Table 1 versus planet mass (Mp, in Mj ), along with published 
1—a error bars for each quantity. For comparison, points for Jupiter and Saturn themselves are also shown. 
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Fig. 2. — Measured planetary radii, Rp (in Rj ), versus central values of the estimated stellar metallicities ([Fe/H]) of the transiting 
planets listed in Table 1, except for WASP-lb and WASP-2b for which metallicity estimates have not yet been published. The names for 
the planets are given in abbreviated form. 
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Fig. 3. — Rp (in Rj ) versus age (in Gyrs) for model planets with masses of 1 Mj (dashed) and 0.3 Mj (solid) for different distances 
[0.02 AU (red), 0.03 AU (yellow), 0.04 AU (green), 0.05 AU (aqua), and 0.06 AU (blue)] from a G2V primary. The models have no cores 
and assume solar metallicites when calculating the opacities. This plot portrays the systematic dependence of irradiated planet radii with 
orbital distance for different masses. See text in 33]for a discussion. 
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Fig. 4. — Solar-opacity-atmosphere/no-core model radii (Rp, in Rj ), at an age of 2.5 Gyr, versus the logarithm base 10 of the stellar flux 
at the planet (Fp), in units of erg cm~^ , for a range of EGP masses from 0.3 Mj to 1.25 Mj . This figure shows both the planet-mass 
and the irradiation-flux dependence of the planet radius, at the average age of stars in the solar neighborhood (~2.5 Gyr). See text in [[3] 
for a discussion. 
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Fig. 5. — The thickness of the radiative zone (A_R), including the transit radius effect, versus mass for coreless models of twelve of 
the transiting planets listed in Table 1. The mean molecular weight (/i) used is that for pure H2/He atmospheres, which is a reasonable 
approximation if the atmospheric heavy-element abundance is not greatly super-solar. Larger /is would translate into smaller Aijs. Since 
age estimates for WASP-lb and WASP-2b are not published, these objects are not included on this plot. The central values of the putative 
ages of the planets are assumed and the calculated thicknesses are given for atmospheric opacities at solar (black) and 10 X solar (red) 
atmospheric values. See text in [[4] for a discussion. 
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10 X Solar, no AR, no core 



OGLE-TR-132b 
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TrES-1 
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Fig. 6. — Rp (in Rj ) versus age (in Gyrs) for a collection of no-core models for the smaller transiting EGPs. They include HD149026b 
(yellow-dashed), HD189744b (green), OGLE-TR-113b (purple-dashed), OGLE-TR-lllb (green-dashed), XO-lb (purple), TrES-1 (gold), 
WASP-2b (blue), and OGLE-TR-132b (red). The top left panel is for solar opacities and does not include the AR term. The top right 
panel is also solar, but does include the AR term. The bottom left panel is for lOxsolar opacities, but does not include the AR term. 
The bottom right panel also assumes lOxsolar opacities, but does include the AR term. This bottom-right panel contains our default 
no-core/no-cloud models. The age of WASP-2b has been arbitrarily set at 2.0±1.0 Gyrs. The barely-perceptible kinks near ~700 Myr in 
the curves for OGLE-TR-132b (red) at the lower left and right and for OGLE-TR-lllb (dashed green) at the lower right are convergence 
glitches in the evolutionary tracks for those models. See discussion in JS] 
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Fig. 7. — Rp (in Rj ) versus age (in Gyrs) for a collection of no-core models for the larger transiting EGPs. They include WASP-lb 
(blue), HATP-lb (aqua), HD209458b (green), TrES-2 (red), OGLE-TR-56b (gold), and OGLE-TR-lOb (yellow). As in Fig. [6] the top left 
panel assumes solar opacities and does not include the AR term. The top right panel is also solar opacities, but does include the AR term. 
The bottom left panel is for 10 X solar atmospheric opacities, but does not include the AR. The bottom right panel also assumes 10 X solar 
opacities, but does include the AR term. This bottom-right panel contains our default no-core/no-cloud models. The age of WASP-lb has 
been arbitrarily set at 2.0±1.0 Gyrs. See [[5]for a discussion. 
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Fig. 8.— Theoretical Rp (in Rj ) versus core mass, Mc (in Earth masses), for OGLE-TR-lOb (yellow), OGLE-TR-56b (gold), HD189733b 
(green), and XO-lb (purple). The lines are for solar, Sxsolar, and lOxsolar atmosphere models and the 3xsolar models are dashed. 
Central values of the estimated stellar ages (Table 2) are assumed. The measured radii of these transiting EGPs, along with l-a error bars 
(vertical), are given. The dots are put arbitrarily at core masses that represent Sxsolar metallicity for the given EGP's measured mass and 
the rightmost extent of the horizontal "error bars" is placed at 3x stellar metallicity masses. If the central value of the estimated stellar 
metallicity is below solar (as for HD189733b), the line is truncated at the dot. Note that to construct the dots the heavy-element fractions 
of the atmosphere and of the envelope/core are here set equal. See text in [|6]for explanations and a discussion. 
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Fig. 9. — Estimated core masses, Mc, (dots, in Earth masses) versus the measured stellar metallicities for twelve of the transiting EGPs 
listed in Table 1. The values are taken from Table 3, where the mean estimated ages of the systems (Table 2) are assumed. For each 
EGP, values for solar (red), 3xsolar (green), and lOxsolar (blue) opacities are given. Since stellar metallicities for WASP-lb and WASP-2b 
are not published, these EGPs are not included on this plot. Note that despite the clustering at low core masses for the large-radius 
exemplars (particularly HD209458b and HAT-Plb) and the low core masses for the moderate-stellar-metallicity EGPs OGLE-TR-56b and 
OGLE-TR-lOb, there is a roughly linear (or, better, monotonic) correlation between metallicity and estimated core mass. The gold points 
indicate the approximate core masses necessary to fit the measured radii when the EGPs in question boast an extra internal power equal 
to a fixed 0.3% of the corresponding hp (Table 4). For these last models, solar-opacity atmospheres, but no irradiation or AR terms, are 
presumed, i.e., these are toy isolated models with cores and internal heat sources. See text in ^and [[7]for discussions. 



22 





Fig. 10. — Top: Theoretical radii of HD209458b (left) and HAT-Plb (right), in units of Rj , versus age (in Gyrs) for different values 
of an hypothesized core power for irradiated atmospheres with solar-metallicity opacities and no solid inner core. The core power lines 
are identified by the fraction of hp, the total stellar power intercepted by the planet. For HD209458b, this is ~7.86xl0~^ Lq and for 
HAT-Plb it is ~5.29xl0~^ (Table 4). The dashed lines are evolutionary trajectories for the respective isolated EGPs with solar- 
metallicity atmospheres without irradiation and without the AR term. The lines identified with the words "Solar, No core" are the 
no-core/irradiated/solar-atmosphere models of Fig. |8] Bottom: Same as for the top, but for isolated atmospheres with solar opacities, 
without irradiation, and with the indicated core powers. Note that the core powers required in this case are much larger than in the 
irradiated case depicted in the top two panels. See JTlfor a discussion. 



